Genetic diversity of vector-borne zoonotic pathogens in companion dogs and cats, Tianjin, China.

Background: Dogs and cats are the hosts of many vector-borne human pathogens that can be transmitted to humans. Given their direct and intimate contact with humans, companion dogs and cats are considered direct sentinels of vector-borne human pathogens. However, limited information is currently available regarding canine and feline zoonotic pathogens in China. This study detected canine and feline vector-borne human pathogens to better understand the potential risk to humans. Methods: Blood samples were collected from 275 domestic companion animals (117 dogs and 158 cats) living in Tianjin city, China, and the presence of DNA from Anaplasma, Babesia, Bartonella, and Rickettsia was detected by semi-nested polymerase chain reaction (PCR). The PCR products of the expected size were sequenced, and these newly generated sequences were subjected to BLASTN, nucleotide identity, and phylogenetic analyses. Results: A total of 24 blood samples tested positive for vector-borne pathogens in companion dogs and cats in Tianjin city, China, with a relatively low positive rate of 8.7%. Specifically, seven human pathogens, including Rickettsia raoultii, Candidatus Rickettsia jingxinensis, Rickettsia sibirica, Rickettsia felis, Babesia venatorum, Bartonella tribocorum, and Bartonella Henselae, were identified. In addition, Anaplasma ovis with zoonotic potential and Candidatus A. cinensis were detected. Conclusion: Our results indicate substantial genetic diversity in the vector-borne human pathogens circulating in companion dogs and cats. Interventions based on "One Health" should be taken to reduce the potential risks of contracting infection from companion dogs and cats in Tianjin, China.

Construction of a gene model related to the prognosis of patients with gastric cancer receiving immunotherapy and exploration of COX7A1 gene function.

BACKGROUND: GC is a highly heterogeneous tumor with different responses to immunotherapy, and the positive response depends on the unique interaction between the tumor and the tumor microenvironment (TME). However, the currently available methods for prognostic prediction are not satisfactory. Therefore, this study aims to construct a novel model that integrates relevant gene sets to predict the clinical efficacy of immunotherapy and the prognosis of GC patients based on machine learning. METHODS: Seven GC datasets were collected from the Gene Expression Omnibus (GEO) database, The Cancer Genome Atlas (TCGA) database and literature sources. Based on the immunotherapy cohort, we first obtained a list of immunotherapy related genes through differential expression analysis. Then, Cox regression analysis was applied to divide these genes with prognostic significancy into protective and risky types. Then, the Single Sample Gene Set Enrichment Analysis (ssGSEA) algorithm was used to score the two categories of gene sets separately, and the scores differences between the two gene sets were used as the basis for constructing the prognostic model. Subsequently, Weighted Correlation Network Analysis (WGCNA) and Cytoscape were applied to further screen the gene sets of the constructed model, and finally COX7A1 was selected for the exploration and prediction of the relationship between the clinical efficacy of immunotherapy for GC. The correlation between COX7A1 and immune cell infiltration, drug sensitivity scoring, and immunohistochemical staining were performed to initially understand the potential role of COX7A1 in the development and progression of GC. Finally, the differential expression of COX7A1 was verified in those GC patients receiving immunotherapy. RESULTS: First, 47 protective genes and 408 risky genes were obtained, and the ssGSEA algorithm was applied for model construction, showing good prognostic discrimination ability. In addition, the patients with high model scores showed higher TMB and MSI levels, and lower tumor heterogeneity scores. Then, it is found that the COX7A1 expressions in GC tissues were significantly lower than those in their corresponding paracancerous tissues. Meanwhile, the patients with high COX7A1 expression showed higher probability of cancer invasion, worse clinical efficacy of immunotherapy, worse overall survival (OS) and worse disease-free survival (DFS). CONCLUSIONS: The ssGSEA score we constructed can serve as a biomarker for GC patients and provide important guidance for individualized treatment. In addition, the COX7A1 gene can accurately distinguish the prognosis of GC patients and predict the clinical efficacy of immunotherapy for GC patients.

Characteristics of Carcinoembryonic antigen related cell adhesion molecules and their relationship to cancer.

Carcinoembryonic antigen related cell adhesion molecules (CEACAMs), such as carcinoembryonic antigen (CEA) and the oncofetal glycoprotein family, are tumor markers. The CEACAMs consist of 12 different human CEACAMs and 5 different murine CEACAMs. The CEACAM family of proteins participates in multiple biological processes that include the immune response, angiogenesis, and cancer. CEACAMs play a significant role in cancer initiation and development. Increasing evidence suggests that family members may be new cancer biomarkers and targets in that CEACEAMs tend to be aberrantly expressed and therefore may have potential diagnostic and therapeutic importance. This review systematically summarizes the biogenesis, biological properties, and functions of CEACAMs, with a focus on their relationship with cancer and potential clinical application. As our knowledge of the relationships among CEACAMs and cancer increases, and as our understanding of the involved molecular mechanisms improves, new therapeutic strategies will evolve for cancer prevention and treatment of cancer patients.

On the sign of the linear magnetoelectric coefficient in Cr2O3.

We establish the sign of the linear magnetoelectric (ME) coefficient,α, in chromia, Cr2O3. Cr2O3is the prototypical linear ME material, in which an electric (magnetic) field induces a linearly proportional magnetization (polarization), and a single magnetic domain can be selected by annealing in combined magnetic (H) and electric (E) fields. Opposite antiferromagnetic (AFM) domains have opposite ME responses, and which AFM domain corresponds to which sign of response has previously been unclear. We use density functional theory (DFT) to calculate the magnetic response of a single AFM domain of Cr2O3to an applied in-plane electric field at zero kelvin. We find that the domain with nearest neighbor magnetic moments oriented away from (towards) each other has a negative (positive) in-plane ME coefficient,α⊥, at zero kelvin. We show that this sign is consistent with all other DFT calculations in the literature that specified the domain orientation, independent of the choice of DFT code or functional, the method used to apply the field, and whether the direct (magnetic field) or inverse (electric field) ME response was calculated. Next, we reanalyze our previously published spherical neutron polarimetry data to determine the AFM domain produced by annealing in combinedEandHfields oriented along the crystallographic symmetry axis at room temperature. We find that the AFM domain with nearest-neighbor magnetic moments oriented away from (towards) each other is produced by annealing in (anti-)parallelEandHfields, corresponding to a positive (negative) axial ME coefficient,α∥, at room temperature. Sinceα⊥at zero kelvin andα∥at room temperature are known to be of opposite sign, our computational and experimental results are consistent.

Field-induced bound-state condensation and spin-nematic phase in SrCu2(BO3)2 revealed by neutron scattering up to 25.9 T.

In quantum magnetic materials, ordered phases induced by an applied magnetic field can be described as the Bose-Einstein condensation (BEC) of magnon excitations. In the strongly frustrated system SrCu2(BO3)2, no clear magnon BEC could be observed, pointing to an alternative mechanism, but the high fields required to probe this physics have remained a barrier to detailed investigation. Here we exploit the first purpose-built high-field neutron scattering facility to measure the spin excitations of SrCu2(BO3)2 up to 25.9 T and use cylinder matrix-product-states (MPS) calculations to reproduce the experimental spectra with high accuracy. Multiple unconventional features point to a condensation of S = 2 bound states into a spin-nematic phase, including the gradients of the one-magnon branches and the persistence of a one-magnon spin gap. This gap reflects a direct analogy with superconductivity, suggesting that the spin-nematic phase in SrCu2(BO3)2 is best understood as a condensate of bosonic Cooper pairs.

Bacteroides thetaiotaomicron ameliorates mouse hepatic steatosis through regulating gut microbial composition, gut-liver folate and unsaturated fatty acids metabolism.

Gut microbiota plays an essential role in the progression of nonalcoholic fatty liver disease (NAFLD), making the gut-liver axis a potential therapeutic strategy. Bacteroides genus, the enriched gut symbionts, has shown promise in treating fatty liver. However, further investigation is needed to identify specific beneficial Bacteroides strains for metabolic disorders in NAFLD and elucidate their underlying mechanisms. In this study, we observed a positive correlation between the abundance of Bacteroides thetaiotaomicron (B. theta) and the alleviation of metabolic syndrome in the early and end stages of NAFLD. Administration of B. theta to HFD-fed mice for 12 weeks reduced body weight and fat accumulation, decreased hyperlipidemia and insulin resistance, and prevented hepatic steatohepatitis and liver injury. Notably, B. theta did not affect these indicators in low-fat diet (LFD)-fed mice and exhibited good safety. Mechanistically, B. theta regulated gut microbial composition, characterized by a decreased Firmicutes/Bacteroidetes ratio in HFD-Fed mice. It also increased gut-liver folate levels and hepatic metabolites, alleviating metabolic dysfunction. Additionally, treatment with B. theta increased the proportion of polyunsaturated fatty acid in the mouse liver, offering a widely reported benefit for NAFLD improvement. In conclusion, this study provides evidence that B. theta ameliorates NAFLD by regulating gut microbial composition, enhancing gut-liver folate and unsaturated fatty acid metabolism, highlighting the therapeutic role of B. theta as a potential probiotic for NAFLD.

Perioperative toripalimab and chemotherapy in locally advanced gastric or gastro-esophageal junction cancer: a randomized phase 2 trial.

Perioperative chemotherapy is the standard treatment for locally advanced gastric or gastro-esophageal junction cancer, and the addition of programmed cell death 1 (PD-1) inhibitor is under investigation. In this randomized, open-label, phase 2 study (NEOSUMMIT-01), patients with resectable gastric or gastro-esophageal junction cancer clinically staged as cT3-4aN + M0 were randomized (1:1) to receive either three preoperative and five postoperative 3-week cycles of SOX/XELOX (chemotherapy group, n = 54) or PD-1 inhibitor toripalimab plus SOX/XELOX, followed by toripalimab monotherapy for up to 6 months (toripalimab plus chemotherapy group, n = 54). The primary endpoint was pathological complete response or near-complete response rate (tumor regression grade (TRG) 0/1). The results showed that patients in the toripalimab plus chemotherapy group achieved a higher proportion of TRG 0/1 than those in the chemotherapy group (44.4% (24 of 54, 95% confidence interval (CI): 30.9%-58.6%) versus 20.4% (11 of 54, 95% CI: 10.6%-33.5%)), and the risk difference of TRG 0/1 between toripalimab plus chemotherapy group and chemotherapy group was 22.7% (95% CI: 5.8%-39.6%; P = 0.009), meeting a prespecified endpoint. In addition, a higher pathological complete response rate (ypT0N0) was observed in the toripalimab plus chemotherapy group (22.2% (12 of 54, 95% CI: 12.0%-35.6%) versus 7.4% (4 of 54, 95% CI: 2.1%-17.9%); P = 0.030), and surgical morbidity (11.8% in the toripalimab plus chemotherapy group versus 13.5% in the chemotherapy group) and mortality (1.9% versus 0%), and treatment-related grade 3-4 adverse events (35.2% versus 29.6%) were comparable between the treatment groups. In conclusion, the addition of toripalimab to chemotherapy significantly increased the proportion of patients achieving TRG 0/1 compared to chemotherapy alone and showed a manageable safety profile. ClinicalTrials.gov registration: NCT04250948 .

Molecular evidence of Rickettsia canadensis in ticks, Hebei, China.

China was affected severely by tick-borne rickettsiosis, and more than 10 Rickettsia species pathogenic to humans have been identified. In recent years, several Rickettsia members, with unknown pathogenicity, firstly identified abroad have been found in China. In this study, parasitic and questing ticks were recovered from two sampling sites in Hebei, China. Specific primers targeting outer membrane protein B (ompB) gene were designed to test the presence of Rickettsia canadensis by nested Polymerase Chain Reaction (PCR). As a result, a total of 428 ticks, including 232 ticks (including 230 Haemaphysalis longicornis and two H. japonica) from Laiyuan County and 196 (H. longicornis) from Luanping County, were collected. Sequencing of PCR products with the expected size and subsequently BLAST showed that 38H. longicornis ticks tested positive for R. canadensis, with an overall positive rate of 8.8%. In addition, 800-bp ompB gene and nearly complete citrate synthase (gltA) gene were recovered from six randomly selected positive samples to better understand their genetic characteristics. Nucleotide identity and phylogenetic analyses showed that R. canadensis presented geographical clustering with evidence that variants identified in the current study presented closer genetic relationship with others identified in Asian than those found in North America. In addition, epidemiological data suggested that H. longicornis may be the competent vector, and more attention should be paid to R. canadensis due to its zoonotic potential. In sum, R. canadensis was confirmed to be present in Hebei Province, China, and its surveillance in ticks should be strengthened due to potential pathogenicity, higher positive rate in ticks and wide distribution of possible vector tick species.

Development and validation of a nomogram for predicting hospitalization longer than 14 days in pediatric patients with ventricular septal defect-a study based on the PIC database.

Background: Ventricular septal defect is a common congenital heart disease. As the disease progresses, the likelihood of lung infection and heart failure increases, leading to prolonged hospital stays and an increased likelihood of complications such as nosocomial infections. We aimed to develop a nomogram for predicting hospital stays over 14 days in pediatric patients with ventricular septal defect and to evaluate the predictive power of the nomogram. We hope that nomogram can provide clinicians with more information to identify high-risk groups as soon as possible and give early treatment to reduce hospital stay and complications. Methods: The population of this study was pediatric patients with ventricular septal defect, and data were obtained from the Pediatric Intensive Care Database. The resulting event was a hospital stay longer than 14 days. Variables with a variance inflation factor (VIF) greater than 5 were excluded. Variables were selected using the least absolute shrinkage and selection operator (Lasso), and the selected variables were incorporated into logistic regression to construct a nomogram. The performance of the nomogram was assessed by using the area under the receiver operating characteristic curve (AUC), Decision Curve Analysis (DCA) and calibration curve. Finally, the importance of variables in the model is calculated based on the XGboost method. Results: A total of 705 patients with ventricular septal defect were included in the study. After screening with VIF and Lasso, the variables finally included in the statistical analysis include: Brain Natriuretic Peptide, bicarbonate, fibrinogen, urea, alanine aminotransferase, blood oxygen saturation, systolic blood pressure, respiratory rate, heart rate. The AUC values of nomogram in the training cohort and validation cohort were 0.812 and 0.736, respectively. The results of the calibration curve and DCA also indicated that the nomogram had good performance and good clinical application value. Conclusion: The nomogram established by BNP, bicarbonate, fibrinogen, urea, alanine aminotransferase, blood oxygen saturation, systolic blood pressure, respiratory rate, heart rate has good predictive performance and clinical applicability. The nomogram can effectively identify specific populations at risk for adverse outcomes.

Identification and Functional Evaluation of Alternative Splice Variants of Dax1 in Mouse Embryonic Stem Cells.

Dax1 (Nr0b1; Dosage-sensitive sex reversal-adrenal hypoplasia congenital on the X-chromosome gene-1) is an important component of the transcription factor network that governs pluripotency in mouse embryonic stem cells (ESCs). Functional evaluation of alternative splice variants of pluripotent transcription factors has shed additional insight on the maintenance of ESC pluripotency and self-renewal. Dax1 splice variants have not been identified and characterized in mouse ESCs. We identified 18 new transcripts of Dax1 with putative protein-coding properties and compared their protein structures with known Dax1 protein (Dax1-472). The expression pattern analysis showed that the novel isoforms were cotranscribed with Dax1-472 in mouse ESCs, but they had transcriptional heterogeneity among single cells and the subcellular localization of the encoded proteins differed. Cell function experiments indicated that Dax1-404 repressed Gata6 transcription and functionally replaced Dax1-472, while Dax1-38 and Dax1-225 partially antagonized Dax1-472 transcriptional repression. This study provided a comprehensive characterization of the Dax1 splice variants in mouse ESCs and suggested complex effects of Dax1 variants in a self-renewal regulatory network.

Tuning magnetoelectricity in a mixed-anisotropy antiferromagnet.

Control of magnetization and electric polarization is attractive in relation to tailoring materials for data storage and devices such as sensors or antennae. In magnetoelectric materials, these degrees of freedom are closely coupled, allowing polarization to be controlled by a magnetic field, and magnetization by an electric field, but the magnitude of the effect remains a challenge in the case of single-phase magnetoelectrics for applications. We demonstrate that the magnetoelectric properties of the mixed-anisotropy antiferromagnet LiNi1-xFexPO4 are profoundly affected by partial substitution of Ni2+ ions with Fe2+ on the transition metal site. This introduces random site-dependent single-ion anisotropy energies and causes a lowering of the magnetic symmetry of the system. In turn, magnetoelectric couplings that are symmetry-forbidden in the parent compounds, LiNiPO4 and LiFePO4, are unlocked and the dominant coupling is enhanced by almost two orders of magnitude. Our results demonstrate the potential of mixed-anisotropy magnets for tuning magnetoelectric properties.

Understanding unconventional magnetic order in a candidate axion insulator by resonant elastic x-ray scattering.

Magnetic topological insulators and semimetals are a class of crystalline solids whose properties are strongly influenced by the coupling between non-trivial electronic topology and magnetic spin configurations. Such materials can host exotic electromagnetic responses. Among these are topological insulators with certain types of antiferromagnetic order which are predicted to realize axion electrodynamics. Here we investigate the highly unusual helimagnetic phases recently reported in EuIn2As2, which has been identified as a candidate for an axion insulator. Using resonant elastic x-ray scattering we show that the two types of magnetic order observed in EuIn2As2 are spatially uniform phases with commensurate chiral magnetic structures, ruling out a possible phase-separation scenario, and we propose that entropy associated with low energy spin fluctuations plays a significant role in driving the phase transition between them. Our results establish that the magnetic order in EuIn2As2 satisfies the symmetry requirements for an axion insulator.

Parallel Genome-Wide CRISPR Screens to Identify State-Dependent Self-Renewal Regulators of Mouse Embryonic Stem Cells.

The pluripotency of embryonic stem cells (ESCs) is more accurately viewed as a continuous developmental process rather than a fixed state. However, the factors that play general or state-specific roles in regulating self-renewal in different pluripotency states remain poorly defined. In this study, parallel genome-wide CRISPR/Cas9 knockout (KO) screens were applied in ESCs cultured in the serum plus LIF (SL) and in the 2i plus LIF (2iL) conditions. The candidate genes were classified into seven groups based on their positive or negative effects on self-renewal, and whether this effect was general or state-specific for ESCs under SL and 2iL culture conditions. We characterized the expression and function of genes in these seven groups. The loss of function of novel pluripotent candidate genes Usp28, Zfp598, and Zfp296 was further evaluated in mouse ESCs. Consistent with our screen, the knockout of Usp28 promotes the proliferation of SL-ESCs and 2iL-ESCs, whereas Zfp598 is indispensable for the self-renewal of ESCs under both culture conditions. The cell phenotypes of Zfp296 KO ESCs under SL and 2iL culture conditions were different. Our work provided a valuable resource for dissecting the molecular regulation of ESC self-renewal in different pluripotency states.

Construction and validation of a prognostic prediction model for gastric cancer using a series of genes related to lactate metabolism.

Background: Gastric cancer (GC) is one of the most common clinical malignant tumors worldwide, with high morbidity and mortality. The commonly used tumor-node-metastasis (TNM) staging and some common biomarkers have a certain value in predicting the prognosis of GC patients, but they gradually fail to meet the clinical demands. Therefore, we aim to construct a prognostic prediction model for GC patients. Methods: A total of 350 cases were included in the STAD (Stomach adenocarcinoma) entire cohort of TCGA (The Cancer Genome Atlas), including the STAD training cohort of TCGA (n = 176) and the STAD testing cohort of TCGA (n = 174). GSE15459 (n = 191), and GSE62254 (n = 300) were for external validation. Results: Through differential expression analysis and univariate Cox regression analysis in the STAD training cohort of TCGA, we screened out five genes among 600 genes related to lactate metabolism for the construction of our prognostic prediction model. The internal and external validations showed the same result, that is, patients with higher risk score were associated with poor prognosis (all p < 0.05), and our model works well without regard of patients' age, gender, tumor grade, clinical stage or TNM stage, which supports the availability, validity and stability of our model. Gene function analysis, tumor-infiltrating immune cells analysis, tumor microenvironment analysis and clinical treatment exploration were performed to improve the practicability of the model, and hope to provide a new basis for more in-depth study of the molecular mechanism for GC and for clinicians to formulate more reasonable and individualized treatment plans. Conclusions: We screened out and used five genes related to lactate metabolism to develop a prognostic prediction model for GC patients. The prediction performance of the model is confirmed by a series of bioinformatics and statistical analysis.

Identification and functional comparison of novel alternatively spliced isoforms of human YAP.

As a key effector of the Hippo pathway, yes-associated protein (YAP) is a major regulator of cell proliferation and apoptosis. In this study, 23 hYAP isoforms were identified in HEK293 cells, with 14 isoforms being reported for the first time. These isoforms were classified into hYAP-a and hYAP-b isoforms based on the variation in exon 1. The two groups of isoforms showed distinctly different subcellular localizations. hYAP-a isoforms could activate TEAD- or P73-mediated transcription, affect the proliferation rate, and enhance the cellular chemosensitivity of HEK293 cells. Moreover, different activation abilities and pro-cytotoxic effects were observed among hYAP-a isoforms. However, hYAP-b isoforms were not found to exert any significant biological effects. Our findings add to the knowledge of YAP gene structure and protein-coding capacity and will help in the elucidation of the function and related molecular mechanisms of the Hippo-YAP signaling pathway.

TIPE2 knockdown exacerbates isoflurane-induced postoperative cognitive impairment in mice by inducing activation of STAT3 and NF-κB signaling pathways.

Objective: Anesthetic exposure causes learning and memory impairment, the mechanisms of which remain unknown. It has been reported that tumor necrosis factor-α-inducer protein 8-like 2 (TIPE2) is a newly discovered immune negative regulator that is essential for maintaining immune homeostasis. This study aimed to examine the role of TIPE2 in isoflurane-induced postoperative cognitive decline (POCD). Methods: An AAV empty vector and AAV shTIPE2 vector for the knockdown of TIPE2 were injected into the dorsal hippocampus of mice. Mice were continuously exposed to 1.5% isoflurane followed by abdominal exploration. Behavioral tests including the open field test and fear conditioning test were performed on the third and fourth day post-operation. Apoptosis was detected by terminal deoxynucleotidyl-transferase-mediated dUTP nick end labeling staining. The kits were used to detect the activity of antioxidant enzymes. Inflammatory cytokine levels were detected by enzyme-linked immunosorbent assay. Signal transducer and activator of transcription 3 (STAT3) and nuclear factor-κB (NF-κB) signaling pathway activities were detected by western blotting. Results: TIPE2 expression increased after isoflurane anesthesia and surgery. TIPE2 deficiency aggravated cognitive impairment in mice and further caused apoptosis and oxidative stress in hippocampal neurons. TIPE2 deficiency induced microglial activation and increased secretion of proinflammatory cytokines. In addition, TIPE2 deficiency promoted STAT3 and NF-κB signaling activation induced by isoflurane anesthesia and after surgery. Conclusion: TIPE2 may play a neuroprotective role in POCD by regulating STAT3 and NF-κB pathways.

The role of KLF transcription factor in the regulation of cancer progression.

Kruppel-like factors (KLFs) are a family of zinc finger transcription factors that have been found to play an essential role in the development of various human tissues, including epithelial, teeth, and nerves. In addition to regulating normal physiological processes, KLFs have been implicated in promoting the onset of several cancers, such as gastric cancer, lung cancer, breast cancer, liver cancer, and colon cancer. To inhibit cancer progression, various existing medicines have been used to modulate the expression of KLFs, and anti-microRNA treatments have also emerged as a potential strategy for many cancers. Investigating the possibility of targeting KLFs in cancer therapy is urgently needed, as the roles of KLFs in cancer have not received enough attention in recent years. This review summarizes the factors that regulate KLF expression and function at both the transcriptional and posttranscriptional levels, which could aid in understanding the mechanisms of KLFs in cancer progression. We hope that this review will contribute to the development of more effective anti-cancer medicines targeting KLFs in the future.

A multi-omics integrative analysis based on CRISPR screens re-defines the pluripotency regulatory network in ESCs.

A comprehensive and precise definition of the pluripotency gene regulatory network (PGRN) is crucial for clarifying the regulatory mechanisms in embryonic stem cells (ESCs). Here, after a CRISPR/Cas9-based functional genomics screen and integrative analysis with other functional genomes, transcriptomes, proteomes and epigenome data, an expanded pluripotency-associated gene set is obtained, and a new PGRN with nine sub-classes is constructed. By integrating the DNA binding, epigenetic modification, chromatin conformation, and RNA expression profiles, the PGRN is resolved to six functionally independent transcriptional modules (CORE, MYC, PAF, PRC, PCGF and TBX). Spatiotemporal transcriptomics reveal activated CORE/MYC/PAF module activity and repressed PRC/PCGF/TBX module activity in both mouse ESCs (mESCs) and pluripotent cells of early embryos. Moreover, this module activity pattern is found to be shared by human ESCs (hESCs) and cancers. Thus, our results provide novel insights into elucidating the molecular basis of ESC pluripotency.